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Abstract. Within the framework of loop quantum cosmology, there exists a semi- 
classical regime where spacetime may be approximated in terms of a continuous 
manifold, but where the standard Friedmann equations of classical Einstein gravity 
receive non-perturbative quantum corrections. An approximate, analytical approach 
to studying cosmic dynamics in this regime is developed for both spatially flat and 
positively-curved isotropic universes sourced by a self-interacting scalar field. In 
the former case, a direct correspondence between the classical and semi-classical 
field equations can be established together with a scale factor duality that directly 
relates different expanding and contracting universes. Some examples of non-singular, 
bouncing cosmologies are presented together with a scaling, power-law solution. 
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1. Introduction 

The old idea pQ of bouncing and oscillating universes has attracted renewed interest in 
recent years. Leading alternatives to inflationary cosmology, such as the pre-big bang 
scenario (for a review, see, e.g., Ref. j2j Ej) and the cyclic/ekpyrotic model |U 151 E] 
developed from brane dynamics in heterotic M-theory IH1 El HHj- both involve a 
contracting universe that undergoes a bounce into an expanding phase. 

However, the physical process that leads to the bounce in these scenarios is unclear. 
Indeed, within the context of classical Einstein gravity, a contracting, spatially flat 
Friedmann-Robertson- Walker (FRW) universe can only undergo a non-singular bounce 
into an expansionary phase if the null energy condition is violated jTTj. Exotic forms 
of 'phantom' matter with negative kinetic energy are often invoked to overcome this 
shortcoming of the classical theory [121 [TBI E! > even though this type of matter leads 
to negative norm states at the quantum level. Quantum gravitational effects can also 
lead to non-singular cosmologies. Models inspired by string/M-theory include those 
where dilatonic loop and higher-order curvature corrections to the tree-level action are 
introduced ^51 QUI Ej- A number of bouncing braneworld scenarios have also been 
developed |TH1 EEE1 113 EU 122! • 

An alternative leading candidate for a non-perturbative and background- 
independent theory of quantum gravity is loop quantum gravity (LQG). (For reviews, 
see, e.g., Ref. [2*3*1 \'24\). This is a canonical quantization of Einstein gravity in terms 
of Ashtekar variables [25! ■ Loop quantum cosmology (LQC) is the restriction of LQG 
to symmetric states. Below a critical value of the scale factor, a.j, the discrete nature 
of spacetime is important and the Hamiltonian constraint is a difference equation for 
the wavefunction of the universe [""o"| . Classical dynamics is recovered when the size 
of the universe exceeds a second critical value, a*. Depending on the quantization 
scheme employed, there is an intermediate 'semi-classical' regime, Oj < a « a„ 
where the dynamics of the universe is determined by coupled, ordinary differential 
equations (ODEs), but where non-perturbative quantum corrections are important. 
It was recently shown that these corrections generically lead to non-singular bouncing 
and oscillatory behaviour [23 12EJ 123 H3 ED] • 

In view of the above developments, the purpose of the present work is to develop an 
approximate analytical framework for studying the dynamics of FRW cosmologies within 
the context of semi-classical LQC when the universe is sourced by a self-interacting 
scalar field. In Section 2, we briefly review the main features of LQC, highlighting 
the quantization ambiguities that lead to important phenomenological consequences. 
Section 3 develops the analytical framework by identifying the scalar field as a natural 
dynamical variable of the system. In the case of spatially flat models, a formal 
correspondence is established between the quantum-corrected Friedmann equations 
and those of classical Einstein gravity. This enables a 'scale factor duality' that 
relates different backgrounds through an inversion of the scale factor to be identified. 
This duality can be extended to spatially-curved cosmologies. Classes of non-singular 
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bouncing cosmologies are found in Section 4, together with a power-law (scaling) 
solution. We conclude with a discussion in Section 5. 



2. Dynamics in Loop Quantum Cosmology 



The physical classical variables in isotropic cosmology are the triad component \p\ = a 2 
and the connection component c = \{k — ja), where 7 ~ 0.274 is the Barbero- 
Immirzi parameter [321 EH] and k = (0, 1) for a (topologically compact) spatially flat 
or positively-curved universe, respectively!. They form a canonical pair with Poisson 
bracket {c, \p\} = 87r7G/3 and the dynamics is determined by imposing the Hamiltonian 
constraint on these variables: 



H 



Ha 



2nGT 



c(c - k) + (l + 7 2 )y 



\p\ +H (t> = 



-a-% + a A V{4>)i 



(1) 
(2) 



where the matter Hamiltonian, 7i^, corresponds to that of a minimally coupled scalar 
field, 0, with a self-interaction potential, V((f>), and canonically conjugate momentum 
Pet, = a?(j). 

The basic variables in LQG are holonomies along curves in space associated with the 
connection and triad. For FRW models, the holonomies are given by hi = exp(cr,) = 
cos(c/2) + 2rjSin(c/2), where = — icrj/2 e SU(2) and <jj are the Pauli matrices. 
Quantization of constraint follows by representing the holonomies as operators and 
promoting the Poisson bracket to a commutator. A primary consequence of LQC 
is that operators for inverse powers of the spatial volume have eigenvalues that are 
bounded and do not diverge at zero volume [HE]- The matter Hamiltonian (J2J) contains 
a divergent factor of a" 3 . This is quantized by starting from the classical identity 
a~ 3 = [3(87T7G/) _1 {c, |p|'}] 3// ( 2-2/ ), where / is a constant, and replacing the connection 
with holonomies j37] : 



4tt7G/ 



^(rjhjihj 1 ^ 21 } 



3/(2-20 



-^tr.^M^V'} 



3/(2-2Z) 



_87T7G/j(j + l)(2j + l) z r^^ J '" L '" J 'J ' (3) 

where the second expression follows because any irreducible SU(2) representation with 
spin j may be chosen when computing the trace jSHEHl- in general, j takes half-integer 
values and j = 1/2 corresponds to the fundamental representation. Eq. (jH)) is a classical 
identity, but inverse powers of the scale factor are not needed if < / < 1. 

I LQG is a canonical quantization and therefore requires a Hamiltonian formulation of the dynamics. 
The negatively-curved FRW cosmology represents the isotropic limit of the Bianchi type V universe 
and, since this is an example of a Bianchi class B model |34|. its evolution can not be described in 
terms of a standard Hamiltonian treatment . We therefore do not consider this model further. 
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Quantization results in a discrete spectrum for the geometrical density operator, 
djj = a~ 3 , that is well-defined and can be approximated in terms of a continuous 
function of the scale factor such that djj(a) = Di(q)a~ 3 , where |3*7| l3*Hj 

DM) = {ji ql ~ l I (/ + 2rl + 1 

- y^*? ((<? + l) m " sgn(g - - IT 1 

and g = a 2 /a 2 , a\ = a 2 j/3, = y/^ylpi and £pi = V GTi is the Planck length. We refer to 
Eq. (J3| as the 'eigenvalue function'. The discrete nature of spacetime becomes important 
below the scale Oj. In the limit a,i < a C o», the eigenvalue function asymptotes to 

D t (a) oc (3/(1 + l)) 3 ^ 2 - 2l \a/a^ 2 - l)/il - l) (5) 

whereas classical behaviour, corresponding to d^i oc a~ 3 , is recovered for a > a*. The 
eigenvalue function is peaked around a ~ a*. 

Hence, the geometrical density differs significantly from its classical limit§ below 
the critical scale a*. Moreover, the freedom in writing the same classical expression 
(JHJ) with different values of (J, I) leads to ambiguities at the quantum level that have 
important cosmological consequences jlH E3 EHl EH EOl H21 E3] ■ For j > 3, there exists 
a semi-classical regime (a^ < a < a*), where spacetime may be viewed as a continuum 
but quantum effects in the density remain important. The dynamics in this regime is 
determined by replacing the a~ 3 term in the matter Hamiltonian (j2J) with dj^{a) 



H = -g^- (a 2 + k 2 ) a + itZ^ + aV = 0, (6) 
where = djfcj). 

In the following Section we develop an approximate analytical approach to FRW 
cosmologies within the context of semi-classical LQC. 



3. Scalar Field Dynamics in Loop Quantum Cosmology 



The Hamiltonian © determines the dynamics via the constraint equation 7i - 
the equations of motion c = — {c, 7i\ and p^ = —{p^.Ti}. These take the form 



H 2 




and 



(7) 



- V 



(8) 



§ In the spatially flat classical model, the Friedmann equations are invariant under a rescaling of the 
scale factor, a — > Ca, for some arbitrary constant C, and it is conventional to normalize such that a = 1 
at the present epoch. The question then arises as to the physical significance of the scale a*. In the 
quantum theory, the primary variable is the volume operator and the scale factor arises by taking the 
cube root. In the quantum theory, therefore, a different normalization for the scale factor is adopted 
and it is the ratio a*/ai that determines the range of the semi-classical regime 00] • The physical scale 
is not specified until has been evaluated and this can only be done at the quantum level. 
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respectively, where H = a/a defines the Hubble parameter and a subscript '0' denotes 
differentiation with respect to the scalar field. 

In general, the dynamics in the semi-classical regime is complicated due to the 
non-trivial form of the eigenvalue function (Jlj). An important point, however, is that 
the transition in going from the semi-classical (a < a,) to classical (a > a*) domains is 
extremely rapid: numerical results indicate that once the eigenvalue function approaches 
unity it does so very quickly and the asymptotic power-law © for the eigenvalue 
function remains a good approximation until a ~ a* |43j . We therefore consider the 
dynamics in the region where Eq. (jSJ) is valid. Such an approximation improves as 
progressively higher energy scales are considered. It then proves convenient to choose 
units such that £p 2 = 8vr and to define the constant = [3/(1 + l)] 3 / (2 - 2 ^al (l ~ 2)/{1 ' l) . 
Eqs. 0-© are then equivalent to 

1 r/) 2 k 2 

ID* a r a z 

where r = 3(2 — — I) > 6 and q = (r — 6)/6. 

We now proceed to consider spatially flat and positively-curved cosmologies in turn. 



3.1. Spatially Flat Cosmologies and Scale Factor Duality 

For spatially flat cosmologies, introducing the new set of variables: 

maps Eqs. (|T0|) and (|TT| into the classical Einstein field equations sourced by a minimally 
coupled, self-interacting scalar field, if): 

3f3 2 = ^ 2 + W (13) 

/5 = -^ 2 , (14) 

where /3 = b/b = —qH. The parameters (6, j3) may be viewed as the rescaled scale factor 
and Hubble parameter, respectively. 

This formal correspondence between classical and (semi-classical) loop quantum 
cosmological scenarios has a number of interesting consequences. It represents an 'ultra- 
violet/infra-red' duality, in the sense that the dynamics of expanding or contracting 
cosmologies in semi-classical LQC can be understood directly by employing methods 
that have been developed previously in the classical scenario. In particular, since q > 0, 
the correspondence relates a classical expanding cosmology with a contracting LQC 
solution, and vice-versa. 
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Moreover, if the scalar field is a monotonically varying function of proper time, Eqs. 
(JT3J) and (JUj) can be transformed into the ' Hamilton- Jacobi' form |3H 1371 

mi EH: 

W(iP) = 3(3 2 - 2(31 ( 15 ) 

= -\i>, (16) 

where a subscript denotes differentiation with respect to the rescaled field, ip, an d all 
variables are viewed as explicit functions of this parameter. It follows from the definition 
of the /3-parameter that 

he* = (17) 

and Eq. (fT7|) may then be integrated to yield the rescaled scale factor in terms of a 
single quadrature: 



b{ip) = exp 



1 H> (3 

2 J d % 



(11 



Expressing the dynamics in this way allows different expanding and contracting 
backgrounds to be directly related through a duality transformation [HU E21- To be 
specific, let us consider a solution to Eqs. (fT5j) and (fTtjj) that is parameterized by the 
functions \j3(i/}), b(ip), W(ip)} and choose as an ansatz for a new solution = b(ip)- 
Substituting this ansatz into Eq. (fT%|) . and noting that the function b(ip) satisfies Eq. 
f!17|) . then implies that the new (rescaled) scale factor and potential are given by 

m = (19) 

q*V^) = 3b 2 (ij)-2bl(ij), (20) 

respectively. On the other hand, substituting an ansatz of the form $(ip) = l/b(ijj) into 
Eq. (fTHjl implies that the dual scale factor is given by b(ip) = 1/ (3(ip) (up to an irrelevant 
constant of proportionality). It then follows from the definitions (|T2*|) that the two new 
solutions are related by a 'scale factor duality' of the form 

a(v) = (21) 

where the corresponding potentials are determined by 

Vty) = V^)b\ii). (22) 

Since this duality relates the scale factor to its inverse, it manifestly pairs an expanding 
cosmology with a contracting universe. It should be emphasized, however, that the 
scalar field tj) does not represent the physical field in the LQC setting. The dependence 
of both the scale factor and potential on the field <fi is determined from the definition 
JEJ such that <p = Dl^q- 1 Jdipa r / 2 (ip). 
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3.2. Positively-Curved Cosmologies 

For positively-curved models, the curvature term in Eq. (fT0|) implies that the 
redefinitions (fT2"|) can not be employed to map the field equations (1TT)|) and (JUJ) onto the 
classical Einstein system. On the other hand, a similar duality to that discussed above 
can be uncovered. Defining an effective energy density 

+ ( 23 ) 
implies that the scalar field equation Q can be written in the form 

a = -3[3ip 2 . (24) 
For a monotonically varying field, this implies that 

^ = -3^ (25) 
and it follows from the definition of the /3-parameter that 

3(3 2 = (26) 
b 

Defining a new variable x = b~ 2 ^ q and substituting Eqs. (|2"3j) and (|2Tj) into Eq. (fTUj) 
then implies that the Friedmann equation can be expressed in the form 

2 

X^ ~ ~X° = -2k q. (27) 

It may be verified that Eq. (}2"Tj) reduces to Eq. (|T7j) in the absence of spatial 
curvature. Since Eq. (|2Tj) is invariant under the simultaneous interchange x(V ; ) cr ('0)) 
we may conclude that a given solution parametrized by [xW, o'W] can be mapped onto 
a dual model such that the new scale factor is given by a(^) = cr 1 / 2 (^) and the effective 
energy density is given by a(ijj) = a 2 {ip). 

In the following Section we present some exact solutions to Eqs. (fTUj) and (fTT|) for 
spatially flat cosmologies. 

4. Inflating and Bouncing Cosmologies 

4-1. Scaling Solution 

The correspondence (|12)) - (fT4"|) between classical cosmic dynamics and loop-corrected 
models implies that a solution to the former can serve as a seed for generating a 
new solution in LQC. The question of whether such solutions are stable to linear 
perturbations can also be determined directly from the stability of the seed solution. 

For example, one solution of particular importance in conventional cosmology is 
the power-law model driven by an exponential potential |H4*j : 

b = t 2 ' x \ V=T ln *> W = ^ r/V ^ (28) 
A A 
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where the kinetic and potential energies of the field scale at the same rate. Employing 
Eq. (|T2"|) implies that the corresponding scaling behaviour in LQC is given by 

r 

2(6 - A 2 ) / r N^'.^/r 

Eq. (|2T?j) represents a scaling solution in the sense that the ratio of the effective kinetic 
and potential energies of the field, <p 2 /(a r V), remains constant. It is interesting that 
when the LQC modifications are significant, the potential is given by a simple power of 
the field. 

Linear perturbations about the classical solution (|28|) decay with eigenvalue m = 
(A 2 — 6)/2 in an expanding universe jSS]. Thus, an expanding (contracting) solution with 
a positive (negative) potential is always stable for A 2 < 6 (A 2 > 6). Similar conclusions 
may therefore be drawn for the LQC solution the stable solutions are represented 

by a contracting cosmology, where A 2 < 6 and the field moves up a positive potential 
away from = 0; or an expanding model, corresponding to a time-reversal of Eq. 
where A 2 > 6, the potential is negative and the field moves towards = 0. This latter 
solution represents a superinflationary cosmology. (Note that although the solution 
becomes singular as t — > 0~, the approximation (jSJ) inevitably breaks down before this 
point is attained). 

4-2. Bouncing Cosmologies 

Eq. implies that a > and it follows that a necessary and sufficient condition for 
a contracting spatially flat universe to undergo a non-singular bounce is that the first 
derivative of the Hubble parameter should remain finite at the point where the sign of 
H changes || . Moreover, the potential of the field must be negative at the bounce and, 
consequently, the scalar field must be a monotonically varying function of proper time. 
The framework outlined in Section 3.1 is therefore ideally suited to determining the 
dynamics of spatially flat, bouncing cosmologies. It proves convenient to express Eqs. 
(jl5|) - (|T5jl in terms of the original variables such that 

V = 3H 2 - ^a r Hl (30) 
q z 

H,= ^t (31) 

2D* or y ' 

The scale factor is then given by 

« r (0) = C+^-/^0-^, (32) 
|| We assume implicitly that the scale factor is positive-definite. 
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where C is a constant of integration, and the time dependence of the scalar field follows 
by integrating Eq. (}3*Tj) : 

t =-*-[* ^ (33) 

Eqs. (I30|) - (}3*3*j) imply that the dynamics of semi-classical LQC is determined 
once the Hubble parameter, H(<ft), has been specified as a function of the scalar field. 
The potential and scale factor follow immediately from Eqs. (130)1 and f!32l) and the 
corresponding time-dependences of the scale factor and scalar field can then be deduced 
after integrating Eq. (|3*3|) and inverting the result. 

In view of this, we now proceed to derive two classes of bouncing models by 
specifying the dependence of the Hubble parameter on the value of the scalar field. 
We begin by considering a Hubble parameter of the form 

H(4>) = Br, (34) 

where B and n are arbitrary constants. Integrating Eq. ()32)1 implies that the scale 
factor is given by 

and Eq. (|3T)J) implies that the potential has the form 

_2CB 2 D,n 2 l2n _ 2 [SB 



q 2 V q 



l + -(n-l) 



A 2/1 



(36) 



For the case where the integration constant C = 0, it can be verified that Eqs. 
and (j33j) correspond to the scaling solution (|2"9"j) discussed above. Eq. (j3l?j) is also 
integrable for the case n = 1 and C ^ and the solution in this case is given by 




(37) 

Q 

where it is understood implicitly that — it/ 2 < t < ir/2. The potential (}3T|) has a simple 
quadratic dependence on the scalar field, and differs from the scaling potential only by 
a negative constant, —2B 2 CD if /q 2 . However, the effect of this constant shift in value 
is significant: the solution now represents a cosmology that contracts to a minimum 
radius at t = 0, at which point the field reaches the maximum of its potential, and then 
bounces into an expanding phase as the field rolls down the other side of the potential. 
The integration constant C determines the shift in the potential and represents the 
minimal value of the scale factor. 

In principle, models with more complicated potentials can be derived in a 
straightforward manner by specifying n ^ 1, although it is not always possible to 
integrate and invert Eq. in these cases. 
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Finally, a second class of model can be found that represents a universe that bounces 
into an asymptotic phase of exponential (de Sitter) expansion. It is given by 



where A is an arbitrary positive constant and the constant of integration in Eq. ()32|) is 
chosen to be C = 1/D*. The potential asymptotes to a positive constant for large \<p\ 
and has a single minimum at <p = 0. The universe bounces into an expanding phase as 
the field passes through this minimum. 

5. Conclusion 

In this paper, we have considered expanding and bouncing universes within the context 
of semi-classical loop quantum cosmology, corresponding to the regime where the 
cosmic dynamics is determined by quantum-modified Friedmann equations. The 
analytic approach developed here compliments previous numerical and qualitative 
studies (23 EHl EHl E0J El ESj and applies in the limit where the eigenvalue function of 
the inverse volume operator has a power-law dependence on the scale factor. 

Such a study is important as it yields insight into the generic nature of these models, 
as well as providing a framework for classifying the different types of possible behaviour 
that may arise. In particular, the modified evolution equations for spatially flat models 
can be recast, after appropriate field redefinitions, into the classical Einstein equations 
for a minimally coupled, self-interacting scalar field. A similar correspondence between 
standard relativistic cosmology and the high-energy limit of the Randall-Sundrum type 
II braneworld [HSj has also been identified recently |HZ1 EHl EHl EDJ and it would be 
interesting to explore possible relationships between these apparently unrelated scenarios 
further. 

Within the LQC framework, such a correspondence is intriguing because it indicates 
that the dynamics at ultra-high energy scales, where the universe is just emerging from 
a truly discrete quantum phase described by a difference equation, can be analyzed in 
terms of low-energy classical models, and vice-versa. By expressing physical parameters 
such as the scale factor and Hubble parameter as functions of the scalar field, we 
identified a form-invariance transformation |H1 E21 EHj that inverts the scale factor 
of the universe. This scale factor duality between two different backgrounds directly 
relates expanding and contracting models and can be extended to the spatially-curved 
models. 

The scale factor duality considered above is different to that arising in quantum 
cosmo logical models derived from lowest-order string effective actions [dT!! IH^T Id^ . 

In this latter class of models, the duality may be viewed as a discrete subgroup of a global 
0(d, d) symmetry of the effective action that arises after toroidal compactification over 
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d spatial dimensions if the dilaton potential has a suitable functional form |fofo|. Thus, 
the duality of string quantum cosmology relates different backgrounds arising from the 
same fundamental Lagrangian. In contrast, however, the scale factor duality of LQC 
relates solutions arising from inequivalent actions, i.e., the transformation (|21j ) -([22 |) acts 
between scalar fields with different self-interaction potentials. Moreover, the scale factor 
duality of LQC becomes apparent when the field equations are expressed as functions 
of an effective scalar field defined by the condition (|12|) . whereas the 0(d, d) symmetry 
of string cosmology is only apparent when the degrees of freedom are viewed explicitly 
as functions of cosmic time A further difference between the two approaches is 

that the field is minimally coupled to Einstein gravity in LQC, whereas the dilaton 
is non-minimally coupled to gravity. In view of the (classical) conformal equivalence 
between scalar-tensor and Einstein theories of gravity, therefore, it would be interesting 
to investigate the LQC approach within the context of scalar-tensor gravity models 
such as the Brans-Dicke theory. Such an approach could yield further insight into the 
similarities and differences between the LQC and string quantum cosmology frameworks. 

Finally, we presented simple models that exhibit scaling and bouncing behaviour 
in the semi-classical regime. One of the principle motivations for developing non- 
singular bouncing cosmologies is that they provide an alternative to the inflationary 
scenario for the formation of large-scale structure in the universe. In principle, density 
perturbations on scales larger than the Hubble radius at the epoch of decoupling can be 
generated during a phase of decelerated contraction jHH EH EH EH] • Since a number of 
issues regarding the propagation of perturbations through the bounce presently remain 
unresolved |7UJ EJ E3 EH] > ^ i s important to develop non-singular bouncing models 
as these can provide a solvable framework for analyzing the evolution of perturbations. 
Moreover, the scaling (power-law) solution has played a central role in the standard 
inflationary scenario and represents one of the few known cosmologies where the density 
perturbation spectrum can be calculated without invoking the slow-roll approximation 
|54[ I74j . It would be interesting to investigate whether the scaling solution derived above 
plays a similar role in determining the evolution of density perturbations generated in the 
semi-classical LQC regime. Ultimately, this could lead to important phenomenological 
constraints on LQC inflationary models. 
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